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A serious difficulty in conventional lattice field theory calculations is the coupling between the chiral and continuum limits. With 
both staggered and Wilson fermions, the chiral limit cannot be realized without first taking the limit of vanishing lattice spacing. 
In this talk, we report on extensive studies of the domain wall formulation of lattice fermions, which avoids this difficulty at the 
expense of requiring that fermion propagators be computed in five dimensions. A variety of results will be described for quenched 
and dynamical simulations at both zero and finite temperature. Conclusions about the benefits of this new method and some new 
physical results will be presented. These results were obtained on the QCDSP machine recently put into operation at Columbia and 
the RIKEN Brookhavcn Research Center. 



1 Introduction 

Important theoretical and algorithmic advances have 
opened a new approach to the numerical study of chi- 
ral symmetry in QCD. The two widely studied lattice 
fermion formalisms, staggered and Wilson fermions, both 
present serious difficulties to the numerical simulation of 
chiral symmetry in lattice QCD. Not only do both of 
these lattice descriptions explicitly break most or all of 
the chiral symmetries present in massless QCD, but they 
also obscure the underlying relationship between topol- 
ogy in the gauge sector and zero modes for the fermions. 
As a result, lattice calculations will typically fail to show 
the full consequences of the Goldstone theorem, current 
algebra, or the 't Hooft solution to the !7a(1) problem, 
without explicit extrapolation to the continuum limit. 
Thus, the physics of chiral symmetry discovered in both 
the 60's and the 70's is corrupted by the usual lattice 
discretization of QCD. 

However, by considering Wilson fermions formulated 
in five dimensions with a large negative mass, tow = 
—mo, Kaplan El showed that it is possible to avoid the 
fermion doubling problem while still achieving complete 
chiral symmetry in the much simpler limit of large lattice 
extent in the new fifth dimension. A further important 
advance was achieved by Narayanan and Neubergera who 
generalized Kaplan's approach and recognized that this 
limit of large fifth dimension not only realized the desired 
continuum chiral symmetry but also created a structure 
of exact fermion zero modes with a close relation to the 
continuum Atiyah-Singer theorem and the physics of the 
axial anomaly. An efficient and well-elaborated version 
of this domain wall method was developedj-and analyzed 
by Shamir and by Shamir and Furmanpla providing a 
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practical and attractive method for large-scale numer- 
ical calculations There has now been considerable ex- 
ploratory numerical work suggesting that this method 
lives up to its promise for both the Schwinger model a 
and quenched QCDL3 and that fermion zero-modes with 
the desired properties are realized HB 

We will describe these issues in somewhat more de- 
tail in Sec 
culations 

pleted computers at Columbia (8,196-nodes 0.4Tflops) 
and the RIKEN Llrookhaven Research Center (12,288- 
nodes, 0.6Tflops)E3 These calculations are designed to i) 
establish the extent to which zero-mode effects can be 
seen in practical simulations with fifth-dimension extent 
L s 10 (Section ||), ii) study the chiral symmetry of 
quenched QCD with careful control over the effects of 
L s , examining both the hadron spectrum and the chiral 
condensate (Section |j) and iii) explore QCD thermody- 
namics using domain wall fermions both as a probe of 
the properties of the pure gauge theory and as a method 
to realize complete flavor symmetry in a lattice study 
of the full QCD phase transition. This final calculation 
represents the best-controlled approach yet available to 
examine the role of the axial anomaly in the full QCD, 
chiral phase transition. 

2 Domain Wall Fermions 

2. 1 Formulation 

Our formulation of domairLjvall fermions follows closely 
that proposed by ShamirtiQ The starting point is the 
standard lattice treatment of the gauge variables as 3 x 3 
special unitary matrices U^n), associated with each link 
in a four-dimensional space-time lattice. Here n locates 
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the site from which the corresponding link extends in 
the positive /j th direction. These link variables enter the 
Wilson gauge action in the usual way. The new features 
of the domain wall treatment appear in the fermion ac- 
tion. Here the fermion field, ift(n, s) is 4-spinor and 3- 
component color vector as usual but now depends on 
both the normal lattice coordinate n and a new fifth co- 
ordinate s, lying in the range < s < L s — 1. The action 
takes the form: 

^DWF= '4>(n,s){(Dy/) 
n 7 s;n' ,s' 

+(D 5 ) s , s ,5 ntn ,}ip(n',s') (1) 
Here D w is the usual Wilson Dirac operator: 

(Dw)n,n' = \ + 7 m M™)m<W-A (2) 

+ (1 - ^)U{n%8 n , nl+li - 2<5„,„,}, 

mo the Wilson mass (but with an unconventional neg- 
ative sign), and D$ an additional fifth-dimension piece, 
diagonal in color: 

(A>)s,y = i{(l + 7 5 )^-i + (l-7 5 )<5.s, S '+i-2 

-m/(l - 7 5 )(5 Si0 (5 s ' : l s -i 

-m / (l + 7 5 )^ Ls _ 1 (5 5 .,o}. (3) 

The combination of the boundaries at s — and s = 
L s — 1 and the negative Wilson mass, — mo, allows un- 
usual massless boundary states to form, localized on the 
s = and s = L s — 1 four-dimensional boundaries of 
the five-dimensional problem. The result is two types of 
states, decaying exponentially as one moves away from 
or L s — 1: massless, right-handed particles bound to 
the right-hand wall and massless, left-handed particles 
bound to the left-hand wall. The mixing between these 
two boundary states, given by the small overlap of the ex- 
ponentially decreasing 5-dimensional wave functions, im- 
plies that these states will actually form a massive four- 
component fermion. They will become literally massless 
only in the limit of infinite separation between the walls, 
L s —> oo, a limit in which this exponentially small over- 
lap vanishes. 

The extra mass term with coefficient m/ in Eq. [| 
explicitly couples the right and left walls and gives these 
boundary states an additional mass proportional to mj. 
This extra term provides an explicit mass which is easily 
adjusted in contrast to the overlap-induced mass whose 
dependence on L s is not so precisely known. We will 
attempt to choose L s and m/ so that the effects of the 
mixing are much smaller than those of m/, so the mass of 
the domain- wall states is explicitly proportional torn/. 



2.2 Chiral Symmetry 

The domain wall fermion formulation described above 
should become a theory of massless fermions in the limit 
L s — > oo. This is easy to see in the free field case and can 
be argued to be true for the interacting theory as wella. 
Heuristically, we observe that if the gauge coupling is suf- 
ficiently small, the relevant gauge configurations will be 
smooth on the lattice scale and interact primarily with 
the low-energy, domain wall states in exactly the fash- 
ion of four-dimensional QCD. Of course, the extent to 
which this is actually true in a particular calculation for 
a specific value of L s , must be studied numerically. 

In the usual Wilson formulation of lattice fermions, 
the axial part of the underlying SU(Nf) x SU(Nf) flavor 
symmetry is strongly broken by the dimension-five, "Wil- 
son term" added to make the doublers heavy. A chiral 
theory is expected to be found in the low energy portion 
of the theory provided the mass-related hopping param- 
eter k is tuned to its critical value. However, in practice, 
the effects of this Wilson term have been so large as to 
obscure the character of the chiral phase transition and 
to serious impede the extraction of a variety of weak ma- 
trix elements whose efficient calculation requires the use 
of chiral symmetry. 

The staggered fermion formulation preserves a sin- 
gle chiral symmetry, making the chiral phase transition 
much easier to simulate. However, this approach con- 
tains lattice artifacts which break the normal vector fla- 
vor symmetries making many quantities more difficult to 
interpret than in the Wilson formulation. 

Both the Wilson and staggered formulations of lat- 
tice fermions are expected to show full, physical flavor 
symmetry only as the continuum limit is approached. In 
contrast, Nf species of domain wall fermions, as formu- 
lated above, support a complete SU(Nf) x SU (Nf) flavor 
symmetry which becomes exact in the possibly much less 
demanding L s — > oo limit. 

2.3 Index Theorem on the Lattice 

An important part of the relativistic quantum physics 
of fermions is the axial anomaly which is required for 
proper understanding of the ir° — > 7 7 decay, the rf mass 
and the order of the QCD phase transition. With both 
Wilson and staggered fermions, the symmetry generated 
by the anomalously non-conserved current is explicitly 
broken by order a or order a 2 terms in the lattice fermion 
action. Thus, without a careful study of the limit of 
vanishing lattice spacing, a — > 0, one is unable to be 
certain whether a given anomalous effect results from a 
lattice artifact or will survive in the continuum limit. 

Happily, this fundamental physical phenomona is 
also represented in a new and potentially more reliable 
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fashion by the domain wall fermion formulation. As de- 
veloped by Narayanan and Neuberger, the domain wall 
formulation supports a variant of the Atiyah-Singer index 
theorem in the L s — » oo limit, even for finite lattice spac- 
ing. In the continuum, this theorem relates the winding 
number of the gauge field background to the number of 
zero modes of the Dirac operator and provides the ex- 
plicit mechanism for gauge field topology to effect the 
physics of quarks. 

To understand the connection between gauge field 
topology and exact zero modes of the Dirac operator, 
for domain wall fermions one begins by viewing the five- 
dimensional fermion path integral in a fixed gauge back- 
ground and for a fixed L s as representing both the deter- 
minant of -Ddwf an d as a quantum mechanical expecta- 
tion value of the transfer matrix T5 for unit translation in 
the fifth dimension in the presence of the s-independent 
background gauge field: 



condensate is easily related to the spectrum of Dirac 
eigenvalues through the Banks-Casher formula: 



Z(L S ) = det{£ D wF} = (0\T 5 Ls \0). 



(4) 



Here, the fermion state |0) is a particular four- 
dimensional fermionic state, determined by the bound- 
ary conditions with half of the available states filled. In 
the limit L s — > 00, the right-most factor of Eq. || projects 
onto the eigenstate of T5 with the largest eigenvalue, |0'). 
Thus, we can conclude that 



lim z ._ 00 det{D D wF} « |<0'|0>| s 



(5) 



The Dirac operator of our five dimensional formulation 
will then have exact zero modes whenever the states |0) 
and |0') have a different number of occupied states (and 
are then, necessarily, orthogonal). This connection of 
exact Dirac zero modes with integers defined from the 
gauge fields is in close analogy to the continuum Atiyah- 
Singer theorem. In the continuum limit, such zero modes 
will correspond to topologically non-trivial gauge config- 
urations. Furthermore, these zero modes should continue 
to exist as, lattice corrections and quantum fluctuation 
are addedla. 

In the next Section, we will examine the question of 
how large L s must be for these expected near-zero modes 
to be recognized. 



p{\, m)d\ 
A 2 + m 2 



(6) 



where p(X) is the density of Dirac eigenvalues correspond- 
ing to the gauge configuration or ensemble of configura- 
tions used to compute ipif) and m is the bare quark mass. 

Eq. ^ implies that the presence of a zero mode, 
p{\) = ZS(X) + ... implies a Z/m divergence in ipip as 
m approaches zero. Such zero-mode effects can be easily 
studied by computing -tpip for a fixed gauge config- 
uration, constructed to be close to a Belavin, Polyakov, 
Schwartz, Tyupkin instanton. Such a configuration will 
have an exact zero mode for the continuum Dirac oper- 
ator. In Figure 0, we show rptp as a function of quark 
mass computed in this fixed, instanton-like background 
using the staggered Dirac operator. Instead of the de- 
sired divergence as mf — > 0, one sees a weak in- 
flection around m/ ~ 10 -2 suggesting that the lattice 
spacing has shifted the zero mode from zero to approxi- 
mately this value. Even this weak signal essentially dis- 
appears if even 10% noise is superimposed on the original 
instanton-like background. 

The behavior for the domain wall Dirac operator in 
this same background, shown in Figure ||, is dramatically 
different. There one sees the expected decrease in tpip as 
nif decreases, until mf ~ 2 10 -2 , at which point a clear 
1 jm f signal is seen, even for the curve corresponding to 
L s as small as 6. For finite L s this 1/mt behavior is 
cut off and -tpip becomes constant asm/ falls below the 
residual mass introduced by direct mixing between the 
two overlapping domain wall states. As can be seen in 
the figure, this effect moves to much smaller m,f as L s in- 
creases. Equally important, the appealing picture shown 
in Figure ^| changes very little as noise is superimposed 
on the instanton-like, background solution. 

Thus, it appears that the domain wall formulation 
will allow us to simulate fermions at fixed lattice spacing 
and not-too-large L s which show both the full SU(2) x 
SU(2) flavor/chiral symmetry and also the important 
connection between lattice topology and fermion zero- 
modes which underlies the physics of the axial anomaly. 



3 Fermion Zero Modes 

In order to investigate the degree to which expected 
fermion zero modes will actually be visible for finite 
L s and non-zero lattice spacing, we have looked for 
zero mode effects in the_background of a discretized, 
instanton-like gauge fieldEZI We do this by evaluating the 
volume-averaged, chiral condensate as a function of the 
explicit quark mass, mf. In the continuum, the chiral 



4 QCD at Zero Temperature 

4-1 Hadronic Spectrum 

As our first test of these ideas we have carried out a 
series of careful, quenched calculations orLSxnall 8 3 x 32 
lattices for a variety of values of L s and m/Ej Using phys- 
ical quark operators whose left- and right-handed com- 
ponents were obtained by evaluating the 5-dimensional 
fermion fields on the s = and s = L s — 1 hyperplanes, 
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Figure 1: The chiral condensate computed on an instanton-like 
background evaluated using staggered fermions. The small inflec- 
tion interrupting the steady decrease of (ipip) with decreasing nif 
is caused by the expected "zero-mode", shifted far from zero by 
this coarse 8 4 lattice. 



we had little difficulty computing hadron masses follow- 
ing the usual methods. In Figure ||, we show the result- 
ing p and nucleon masses, extrapolated to m/ = 0, as a 
function of Ls- Although we have examined quite large 
domain wall separations, up to L s — 48, it is clear from 
the figure, that we could have extracted accurate results 
from the smaller L s = 10 or 16 calculations. 

The corresponding L s behavior of the pion mass is 
shown in Figure ||. In contrast to the case of the nucleon 
and the p, we see significant dependence on L s , even 
at the largest values of L s . This contrast can be easily 
explained if we hypothesize that the largest effect of finite 
L s is to give an additional independent contribution 
to the quark mass. Such an effect will be much more 
visible for the very light pion than for the heavier p or 
nucleon. Using the observed to/ sensitivity (not shown 
here) of m p and m n , one can conclude that the 100% 
change in rn^\ m .—o seen as L s increases from 16 to 48, 
would amount to only a 4% effect on m p \ m .—o, an effect 
hidden by our errors. Of greater potential interest is the 
failure of m n to vanish as L s — > oo. While this may be a 
new feature of the quenched approximation made visible 
by the domain wall formalism, it is very likely a simple, 
finite- volume effect. This non- vanishing contribution to 
to w is an appreciable fraction of to^ only for very small 
quark masses, To/ ~ 0.02. However, for such light quarks 
the product to/ (ipip) V = 0.02 • 0.0019 • 8 3 x 32 = 0.7, a 
clear warning that the Goldstone phenomena will begin 
to be influenced by our finite 8 3 x 32, space-time volume. 
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Figure 2: The chiral condensate computed on the instanton-like 
background of Figure 1, now evaluated using domain wall fermions. 
Results for four different choices of L B are shown: L B = 4 (circles), 
6 (squares), 8 (crosses) and 10 (diamonds). Now the expected, 
zero-mode induced 1/my behavior can be easily seen for L s > 4 
and rrif < 10~ 2 . 




Figure 3: The p and nucleon masses, extrapolated to mj = 0, 
plotted as a function of L s , the lattice extent in the fifth dimension. 
These masses were computed with domain wall height mo = 1.65 
from quenched, (3 = 5.85, 8 3 X 32 configurations. 
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Figure 4: The tt mass, extrapolated to zero quark mass, plotted as a 
function of L s , again for quenched, j3 = 5.85, 8 3 X 32 configurations. 
The failure of to approach zero as L s — > oo, as required by the 
Goldstone theorem, is most likely a simple, finite-volume effect. 



4-2 Chiral Condensate 

Lacking the fermion determinant, the quenched approx- 
imation can potentially yield gauge configurations with 
small or vanishing Dirac eigenvalues. As was discussed 
in Section such configurations could lead to an un- 
physical, 1/m divergence in the chiral condensate. How- 
ever, such behavior could be easily suppressed by lat- 
tice artifacts for staggered fermions and may be reduced 
when "exceptional" configurations are discarded for Wil- 
son fermions. 

We have explicitly searched for such effects in a 
/3 = [5,85, quenched simulation and found them quite 
easily!?! Our results for a 8 3 x 32 lattice are shown in 
Figure |. One sees a clearly divergence for small 

to/. The 1/m/ coefficient is 3.8(3) 10 -6 . If this effect 
comes only from exact zero modes, whose number in- 
creases as V 1 / 2 , we should expect this behavior to be 
less pronounced on larger volumes. In fact, a similar, 
16 3 x 32 calculation also shows this 1/m/ behavior, but 
with a much smaller coefficient, 0.6(1) 10 -6 . Thus, while 
one should worry that conventional quenched chiral ex- 
trapolations do not allow for such behavior, this new 
effect may not be important for quenched calculations 
performed on large physical volumes. 

5 QCD near T c 

Because of the importance of chiral symmetry for the 
QCD phase transition, one of the most interesting ap- 
plications of domain wall fermions may be to studies of 
QCD thermodynamics. 
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Figure 5: The chiral condensate is plotted as a function of quark 
mass from the same quenched, /3 = 5.85 8 3 X 32 calculation. The 
clear behavior is a new, unphysical feature of the quenched 

approximation, resulting from the absence of the fermion determi- 
nant. 



5.1 The quenched Chiral Condensate 

We will first continue the discussioiL. above of possible 
quenched, 1/m/ behavior for T > T C \B Here, one expects 
to see such singular behavior and a 1/m/ term that may 
be non-zero, even in the limit of infinite volumecJ While 
earlier staggered fermion calculations have not seen this 
effcctpSEj it is very visible in the /3 = 5.71 results shown 
in Figure ^. A companion calculation on a larger 32 3 x 4 
lattice shows 1/m/ behavior with the same coefficient. 
Thus, by using a fermion formulation more sensitive to 
the effects of zero modes, we have discovered that the 
quenched chiral condensate shows quite striking behav- 
ior. Rather than vanishing for small to/ in the deconfincd 
region, as would be expected if quenched chiral symme- 
try were restored above /3 C , the chiral condensate diverges 
as 1/m/. In addition, as can be seen in Figure |[ there 
is a non-zero, m/ = intercept even if the 1/m/ term 
is removed. Clearly the old picture of quenched chiral 
symmetry restoration above (3 C is seriously incomplete. 

5.2 SU(2) x SU(2) Chiral Transition 

Now let us examine the physical, 2-fiaYor QCD phase 
transition using domain wall fermionsE-3 We have done 
a series of calculations on 8 3 x 4 lattices with a vari- 
ety of values of /3, L s and domain wall heights. In Fig- 
ure 0, we show the Wilson line expectation value and the 
chiral condensate computed with a domain wall height 
mo = 1.90, L s — 12 and quark mass m/ = 0.1. One sees 
the expected cross-over behavior as /3 increases, suggest- 
ing that a full QCD, thermodynamic calculation using 
domain wall fermions is may well be possible. We have 
seen very similar behavior for mo varying between 1.65 
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Figure 6: The chiral condensate is plotted as a function of quark 
mass, computed using pure gauge configurations on a 16 3 x4 lattice, 
just above the deconfining phase transition for f} = 5.71. There is 
an new 1/m f divergence visible for small mass as well as a non-zero 
constant piece in {jyiji. These appear to be new, volume-independent 
terms, first visible when domain wall fermions are used. 
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and 2.15, suggesting that we are seeing stable, 2-fiavor 
behavior for to = 1.9, although more study of this ques- 
tion is warranted. 

Of special interest is the degree of chiral symmetry 
present on either side of this transition region. This is ex- 
amined in Figure |§| where the limit to/ — > is displayed 
using a series of full QCD calculation with quark masses 
m f = 0.06, 0.10, 0.14 and 0.18. The domain wall formu- 
lation appears to have allowed a calculation, for the first 
time, in which a transition is seen between a phase with 
spontaneously broken chiral symmetry (/3 = 5.20) and 
one with what should be full, restored, 577(2) x SU(2) 
symmetry (/3 = 5.45). Calculations are now underway 
to study this system further, examining the properties 
of the QCD phase transition in greater detail for larger 
spatial volumes. 

5.3 Anomalous Symmetry Breaking 

As a final topic, we attempt to exploit the sensitivity of 
the domain wall fermion formulation to topological ef- 
fects, and examine the degree of anomalous-symmetry 
breaking slightly above the phase transitiontfl The cur- 
rent results of this study are shown in Figure ||. We 
plot the difference of two screening masses, computed at 
[3 = 5.40, just above the transition region, evaluated on a 
16 3 x 4 lattice, using L s — 16 and too = 1.9. Here we are 
subtracting the screening mass mg, extracted from the 
exponential damping of correlation functions computed 
from the operator (t/jftp)^), and computing using 
the operator (i^fjf^ 5 t/j)(x). The 2x2 matrices, r, are the 
usual Pauli, isospin generators. 



Figure 7: The chiral condensate and Wilson line are plotted as a 
function of (3. These values were obtained in a full QCD, Nf = 2, 
rrif =0.1 calculation on an 8 3 X 4 lattice. Behavior characteristic 
of a phase transition with j3 c ~ 5.35 is seen. 
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Figure 8: The chiral limit of the (ipijj) in full, 2-flavor QCD for 
L a = 16 and mo = 1.9 Shown are points below the transition, 
/9 = 5.20, 8 3 X 4(squares) and points above the transition, /3 = 5.45, 
8 3 X 4(diamonds) and 16 3 X 4(crosses). 
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Figure 9: The anomalous difference of the screening masses for the 
7r and 8 are plotted as a function of quark mass. These calculations 
were performed with two flavors of dynamical quarks on a 16 3 X 4 
lattice, at /3 = 5.40, just of above the transition, using L a = 16 and 
mo = 1.9 



Because the S and tt correlation functions are related 
by the anomalous chiral symmetry, any difference be- 
tween these two screening masses in the limit rrif — > 0, is 
an unambiguous measure of anomalous chiral symmetry 
breaking. Since these screening masses are both about 1 
in lattice units, Figure || shows a small (« 5%), but pos- 
sibly significant effect. The present data, shown in the 
figure, give an intercept of 0. 067(1. 8)when fitted to the 
expected quadratic m f dependence. While we cannot be 
certain we have seen the sought-after evidence of anoma- 
lous symmetry breaking (the sort required to predict a 
second-order QCD phase transition), we can assert with 
confidence that the effect is quite small. 

6 Conclusion 

We have demonstrated that the domain wall formulation 
of lattice fermions can be used to study both zero tem- 
perature, quenched hadron spectroscopy and the QCD 
phase transition for both zero and two flavors. This 
new method appears to realize accurate chiral symme- 
try at fixed, non-zero lattice spacing, and does not re- 
quire unreasonably large domain wall separation in the 
fifth dimension. This method offers a very promising new 
approach to the study of chiral symmetry using lattice 
methods. 
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